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^ ' Abstract 



We consider some models in the Kardar-Parisi-Zhang universality 
class, namely the polynuclear growth model and the totally/partially 
asymmetric simple exclusion process. For these models, in the limit of 
large time t, universality of fluctuations has been previously obtained. 
In this paper we consider the convergence to the limiting distributions 
and determine the (non- universal) first order corrections, which turn 
^ ' out to be a non-random shift of order t~^^'^ (of order 1 in microscopic 

units). Subtracting this deterministic correction, the convergence is 
then of order We also determine the strength of asymmetry in 

. the exclusion process for which the shift is zero. Finally, we discuss to 

I what extend the discreteness of the model has an effect on the fitting 

' functions. 

(N 

g ; 1 Introduction and results 

T— I ■ 

A growth model is a stochastic evolution for a height function ht{x), x space, 
t time. In this paper we consider some models in the Kardar-Parisi-Zhang 
^ ■ (KPZ) universality class |2T] in 1+1 dimensions, that are irreversible and have 

^ ■ local growth rules. Moreover, there is a smoothing mechanism ensuring the 

existence of a deterministic limit shape /ima(0 •= ^^^t^oot~^ht{^t), around 
which fluctuations are expected to have some degree of universality. More 
precisely, the height function under the scaling 
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should be a well-defined stochastic process as t — ?■ oo, which means that 
fluctuations are of order t^^'^ and the spatial correlation length scales as t^/^ 
(below we will focus on the one-point distribution only). The limit process 
still depends on initial conditions and is expected to be universal inside sub- 
classes of the KPZ class (e.g., the processes for flat or curved limit shape are 
different). For more details, see the recent reviews [THIIIS] and the references 
therein. 

Theory. On the theoretical side, some solvable models in the KPZ class 
have been analyzed in great detail. Two of the best studied models are 
the polynuclear growth (PNG) model and the (totally/partially) asymmet- 
ric simple exclusion process (TASEP/PASEP), see below for a concise def- 
inition. For example, it was shown [UlT^lEHlllI] that if the initial condi- 
tions generate a curved limit shape, then the limiting distribution function 
of Ft = P{hl'''"'{0) < s) as t — )■ oo is given by the GUE Tracy- Widom dis- 
tribution function, -Fque, occurring flrst in random matrix theory in |37] . 
Recently, the flrst model in the KPZ class, namely the KPZ equation itself, 
it was solved and the -Fque was obtained P|IST | [52 |IMP . 

This solution of the KPZ equation includes the distribution function for 
ht at any time t, not only in the t oo limit. A further aspect noticed 
in [31] is that the correction of Ft from -Fque is of order t~^^^. This means 
that on the original scale, the difference between the height function ht{^t) 
and thraa{^) is of order 1. 

Experiments. Until recently, on the experimental side there were only few 
experiments giving the fluctuations exponent 1/3 [211113]. Besides the diffi- 
culties of having good statistics, one of the main issues in the experimental 
set-up is to have really a local dynamic, and the centering in (11. ip has to 
be obtained experimentally from the measured asymptotic growth velocity. 
In any case, experimental data were not good enough to have more detailed 
information on the scaling exponents, until the recent amazing experiments 
carried out by Takeuchi (see [35] and [36]). Using nematic liquid crystals 
they were able to get accurate statistics that confirmed the fiuctuation and 
correlation exponents, but also the limiting distribution functions and the 
covariance of the processed previously obtained in solvable models. 

As it was the case for the solution of the KPZ equation cited above, also in 
these experiments one could see that the fit between the density of the Tracy- 

^ These works uses the approach of Bertini and Giacomin for the weakly asymmet- 
ric simple exclusion process (WASEP) [6]; a replica approach is in [T2j[T4j[28]; see the 
review [33] for details. 

^The processes are: (a) for non-random initial condition and flat limit shape: the Airyi 
process [Tni[3D], (b) for curved hmit shape: the Airy2 [13 process. 
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Widom distributions and the measurements is good even at relatively small 
time t, but finite size corrections are still visible. They measured also the 
decay of the mean, variance, skewness and kurtosis. In the scaled variables, 
the mean has been seen to decay as t~^^^, while the others as Thus, in 

the unrescaled variables, the mean has a shift of order 1. 

Before explaining the results, we briefly describe the models analyzed in 
this papeiEl. 

(a) PASEP and TASEP. The partially asymmetric simple exclusion pro- 
cess (PASEP) on Z in continuous time is an interacting particles system. At 
any time instant t at most one particle can occupy a site in Z and particles 
try to jump to the right neighboring site with rate p and to the left neighbor- 
ing site with rate q = 1 — p. The jumps are made only if the arrival sites are 
free. This dynamics does not change the order of particles. We label them 
from right to left so that Xk{t) denotes the position of particle with label 
k at time t and Xk{t) > Xk+iit) for all k and t. The two initial conditions 
analyzed in this paper are step initial condition, Xk{0) = —k, k = 1,2, . . ., 
and alternating initial condition, Xk{0) = —2k, A; G Z. When g = (and 
p = 1) the particles can only jump to the right so that the model is totally 
asymmetric and is called TASEP. 

(b) PNG model. The polynuclear growth model describes an interface 
given by a height function x ^ ht{x) G Z, x G M, which is a step function 
with up- and down-steps on the integers and constant in between. The PNG 
dynamics is the following. The up-steps move to the left with unit speed and 
the down-steps to the right with unit speed. On top of it, there are random 
nucleation events, i.e., creation of a pair of up- and down-steps, which then 
follow the deterministic spread to the left/right. In this paper we consider 
ho{x) = for all a; G M and nucleations as follows: for flat PNG, nucleations 
form a Poisson point process in M x with intensity 2, while for PNG 
droplet the nucleations are further restricted to {{x,t), \x\ < t} (see also the 
review [T7] for more details and illustrations). 

Results. A difference between the shift of the mean in the solution of the 
KPZ equation and in the liquid crystal experiment is that they have opposite 
signs. The latter has the same sign as for the TASEP (totally asymmetric) 
as we will show. Since the solution of the KPZ equation can be obtained 
starting from the WASEP (weakly asymmetric), there will be a value of the 
asymmetry for which the mean has no shift (up to (9(t"^/^)). A preliminary 



■^Java animations of these models can be found for the PNG model at 
|htt p : //www-wt . iam.uni-bonn .de/~f errari/ animations/AnimationRSK.h tinl| and at 
Ihttp : //www-wt . iam.uni-bonn .de/~ferrari/animations/ContinuousTASEP .html for 
TASEP. 
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Monte-Carlo simulation indicates that this happens for the PASEP height 
function at the origin for p = c:^ 0.78 [HT] . 



Result 1. In Corollary \4 -31 we show that pc is the solution of 
°° (1-vY 1 

E / n M = o ^ pe = 0.78227 87862... (1.2) 

We first determine an analytic formula for the shift of the distribution of 
a tagged particle (see Proposition 14. ip . The shift turns out to be a function 
of the macroscopic particle number. However, when we switch back to the 
height function representation, the shift becomes again independent of the 
macroscopic position (the ^ in (11. ip ). 

The other results concern the first order correction to the limiting distri- 
bution function and density. Let us illustrate it for the PNG droplet. The 
other cases are analogue, but instead of the GUE one has for example the 
GOE Tracy-Widom distribution. The height function is an integer-valued 
function, i.e. ht{0) G Z, and since we always look at it it at a single point, 
X = 0, we drop the space dependence in our notation and write ht = ht{0). 
It is known that for some explicit constants ci,C2, the rescaled variable 

r ht - cit V . , ^ ,^ 

ht,resc ■= ^ ^^/3 ^ C aS t -> OO, (1.3) 

where C is a random variable with (GUE) Tracy-Widom distribution -Fgue, 

Ft{s) := P(/if,resc < s) ^ Fgve{s) as t ^ OO. (1.4) 

Remark that Ft is piecewise constant over intervals of length 5t := l/(c2t^^^). 
One expects that 

ht,resc = C + vSt + 0{6^) on It := {Z - cit)6t. (1.5) 

where t] is another random variable a priori not independent from (. With 
(11.51) we mean that 

nht,rosc <s) = FiC + v5t + 0{5l) < s) (1.6) 

for s G It- 

But what is the nature of r;? The surprising result is that for the models 
we consider, 77 is a deterministic constant and therefore independent from ( 
(see Section E] for PNG and TASEP, Section 1 for PASEP). This implies the 
following. 
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Result 2. Let us denote by 5t := ^t ^^'^ the discrete lattice width where 
ht,resc lives. There exists a constant rj such that 

Ft{s) = P(Vesc <s)= Fgue(s -vSt) + 0{6l) (1.7) 

for all s ^ It = {1j — Cit)6t. 

For PNG and TASEP the shift does not depend on the chosen macroscopic 
position, but this property is not generic and might depend on the chosen 
observable too, as shown by the result on PASEP. (This non-universality of 
rj is quite intuitive, since 77 is a correction term on the microscopic scale, thus 
model-dependent.) Consequently, by shifting the height function ht by the 
constant rj as in (11. 7p . the convergence of the distribution function to Fgue 
is of order 0{t~'^^^). If r] was not independent from then one would have 
a convergence only of order 0{t~^^^) instead. 

In the domain of random matrices, similar results have been obtained for 
the Gaussian and Laguerre Unitary Ensembles [13l|22]. However, in those 
cases the analyzed random variable were continuous. This differs from the 
random variables of the models considered here, since before rescaling they 
live on Z, while after the rescaling (II. ip they still live on a discrete lattice of 
width 6t. The discreteness becomes irrelevant for the universal statements, 
but at first order it can not be neglected for the fit with the limiting dis- 
tribution function and density. Indeed, the shift needed to have a fit with 
accuracy of order 0{t~'^^^) is not the same for the density as for distribution 
function. In order to see this feature, consider the slightly modified scaling 
of the height function 

ht,resc ■= {ht - cit - a)6t (1.8) 
where a G M is a given constant. Further, set 

Ft{s) := P(Vesc < s), (1.9) 

and 

Ptis):=6r\Ftis)-Ftis-6t)). (1.10) 

Result 3. With the choice a := rj + ^ we have 

Ptis) = F^^^is) + Oi6^). (1.11) 

for all s ^ It := ~ Cit — a)6t. 

These results are discussed in Section [2] and used for the fits of the sim- 
ulations of TASEP in Section [31 
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Remark 1.1. Result [3] is generic and does not depend on the fact that our 
distributions is expressed by a Fredholm determinant, but is a consequence 
of the 0{6^) error for the centered discrete derivative. 

Remark 1.2. With the scahng f ll.Sp . Result [2] writes 

Ft{s) = Fgve{s + '^6t) + 0(6^), s G It, (1.12) 
while the scaling (11. 3p yields 

Ftis) = Fgve{s + - a6t) + 0(6^), s e 1, (1.13) 

and 

Pt{s) = Fl,^^{s-a5t) + 0{5^), sElu (1-14) 
where pt{s) := 6T\Ft{s) - Ft{s - 6t)). 

In view of these results, we carried out a simulation for TASEP with time 
t = 1000. As observable we used a tagged particle. The dots in Figure H] 
represent s t— pt{s) for s G It and a = "Z] + |, which is well approximated 
by the solid hne s i— )■ Fq^-^^s) as predicted by Result [31 As comparison, the 
dashed line is the unshifted density s i— )■ Fqjj^{s — a6t) (see (11.141) ). i.e., the 
fit obtained with a = 0. 

The same applies to the distribution function. The dots in Figure [3] 
are the plot of s H- Ft{s) for s E It and a = r] + ^. The dashed line is the 
predicted limiting distribution function with scaling (II. 3p . be s H- -Fgue(s) = 
-^gue('S + ^St — aSt) (see (I1.13P ). The fit suggested by Result [2] is the solid 
line, s H- -Fgue('S + ^St), which indeed is a better fit. 

In the same way we fit Figures [2] and [T] with the difference that the limiting 
distribution function is s t-j- Fgoe(2s). 

Finally, the shift used in ht^iesc is the same needed to have a convergence 
of the moments, and consequently of the variance, skewness, kurtosis of order 
(9(t~^/^). The following result is discussed in Section [2l2l 

Result 4. We have 

mZesc) = nn + 0(5^) (i.is) 

for all m G N. 

Remark that if rj was not independent from the convergence of the 
variance, skewness, and kurtosis would still be of order 0{t~^^^). 
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2 Strategy and effects of the discreteness 

In this section we present the strategy used to get the results. We discuss 
the effects of the intrinsic discreteness of the models on the fitting functions 
and on the moments, since it is relevant at first order. Finally, we explain 
how to fit data coming from an experiment. 

2.1 On the fitting functions (Results [2] and [3]) 

For the PNG model and the TASEP, the strategy of getting (12. 8p is the 
following. In these cases, the distribution function of ht can be expressed as 
a (discrete) Fredholm determinant with kernel K^, 

P(/it <x) = det(l-i^t),2({,+i,,+2,...}), xeZ. (2.1) 

For some constant a G M, the rescaled random variable 

ht,resc ■■= {ht - cit - a)6t with 6t = c^H'^^^ (2.2) 

lives on Jf := (Z — cit — a)6t. According to the scaling in (12.21) . we define the 
rescaled kernel Kt^rcsc as 

-K't,rcsc(si, S2) ■■= S^'^Kticit + a + cit + a + 826^'^) (2.3) 

so that the distribution function Ft defined by 

Ft{s) ■.= P{ht<cit + s6^^ + a), s G M, (2.4) 

can be written as a Fredholm determinant on £"^{[3 + 6t, 00) (1 ItjStiy) with u 
the point measure on It, 

Ft{s) = det(l — -R't,rcsc)f2([s+5t,oo)n/t,<5ti/) 

= ^ ^ ^ Stdet{Kt,,^sc{xi,Xj))i<ij<n. 

n=0 Xl,...,Xn£Jt 

^Mathematically, we get a weaker result, but to illustrate what really happens let us 
assume that one has (|2.8p . What is missing are explicit bounds on the decay of the 
kernels, which can be obtained by standard asymptotic analysis; the ingredients like the 
steep descent paths are all already contained in previous papers. For TASEP wc illustrate 
the results with a simulation for time t = 1000. 
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Note that Ft and the Fredholm determinant in (12.51) are piecewise constant 
functions, with jumps for values of s in the lattice It- The next step is to 
show that for Si, S2 G It and a well-chosen a G M, 

i^f,resc(Sl, S2) = K{SI, S2) + 6tK^sym{Sl, S2) + 0{S^) (2.6) 

where is a symmetric and -ft'asym an antisymmetric kernel. Then, it follows 
thaJl 

Ftis) = det(l — -ft'j,rcsc)£2([s+<5t,oo)n/t,<5ii/) 

= det(l - K)i2(^[s_^St,oo)nit,Sti') (2.7) 
X (1 + 5iTr((l - XsKxs)-')XsK^ymXs) + 0{6^)) 

where Xs is the projection onto [s + 6t, oo) fl It- The operator under the trace 
is antisymmetric, therefore its trace is zero and 

Ft{s) = det(l - i^)^2([,+5,,^)n/„5,.) (1 + Oi5f)) , seR. (2.8) 
If we denote by 

F{s):=det{l-K)L2(^^s,oc)) (2.9) 

the limiting distribution of Ft{s) taken as a Fredholm on L'^{{s, oo)), then by 
Lemma 12.11 below, we get 

Ft{s) = F{s + l5t) + 0{5^), (2.10) 

for s E It, and by the argument below (that gives Result [3]), one finally 
obtains 

Ms) = F'{s) + Oi6^). (2.11) 

In Section [3] we derive ([2J]) for the PNG model and the TASEP. 

Let us explain how to get Result [3] without the need of Fredholm deter- 
minant representations. Assume that there exists a constant 7 such that 

Ft{s) = F{s + -f6t) + 5lQ{s) + 0{5l), s G h, (2.12) 

with F E and Q E C^. Then using Taylor expansion we readily obtain 

Pt{s) = F'{s) + |(f - (1 - l)')F"{s) + 0{6^), s G (2.13) 

Therefore, if 7 = 1/2, then pt{s) - F'{s) = 0{6t) for s G It, while the 
approximation would be only of order 6t if 7 7^ 1/2. 

In our case, see Corollary 12. 2^ we have 7 = 1/2 which is a consequence of 
the following lemma. 

^On a rigorous level, one needs to verify that (a) the 0(5^) in (|2.6p is an operator with 
with 1-norm of order 0(6^) and (b) (1 — XsKxs)~^XsKasyniXs is trace-class. 
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Lemma 2.1. Assume that the kernel K satisfied 

max{|i^(xi,X2)|,|Wxi,X2)|,|9,9,i^(xi,X2)|} <Ce-'=("i+"^) (2.14) 

for some constants C,c> 0, for all Xi,X2 G (s, oo) and i,j G {1,2}. Let 6t 
be as above the lattice width. Thei^ 

|det(l - K)L2(is+5t/2,oc)) - det(l - i^)^2([.+5,,oo)n/t,5t^) | = OiS^e"''). (2.15) 
Proof. Let us set Jt := fl [s + 6t, oo). Then, we have 

°° f— 1)" 

det(l-i^),2(j,,,,,) = ^^-^ J2 5'^detiKixi,x,)),<ij<n (2.16) 

n=0 Xi,...,Xn£Jt 

and 

°° (—1)" /" 

det(l - K)L2(^(^s^St/2,oo)) = ^ j— / d"x det(i^(xi, Xj))i<ij<n. 

(2.17) 



(s+St/2,oor' 



Equation fl2.15p then follows from Lemma lA.ll with 

/(xi, ...,Xn) := det{K{s + 6t + Xi,s + 6t + Xj))i<ij<n (2.18) 
together with Lemma [A. 2 [ □ 
A straightforward corollary is the following. 

Corollary 2.2. Assume (12.81) and (12. lip to hold. Then, for large t, we have 
(remember that 6t = C2^t~^^^) 

Ft{s) = Fis + ISt) + 0(6^) (2.19) 

for s E It- 

Remark 2.3. An equivalent way would be to consider the scaling (12. 3p 
without the shift by a, i.e., 

Kt,resc{Si, S2) ■■= 6^r^Kt{Cit + Cit + 52^^"^) = i^t,resc(Sl - a6t, S2 - a6t) 

(2.20) 



^With di we mean the derivative with respect to the ith entry of the function. The 
assumption (j2.14p holds for the Airy kernels, see Lemma [A. 31 

''For the Airy kernels it is easy to improve ©(t^^Z-^e"") to 0(t-2/3g- max{s,o}^^ However, 
getting a rigorous good bound for the error as s — >■ —00 is a much more difficult task (this 
would be needed for a rigorous proof of the convergence of the moments). 
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Then, instead of fl2.6p we would have obtained 

-f^f,rosc(si, S2) =K{si, S2) - a6t{diK{si, S2) + d2K{si, S2)) 

+ 6tK^y,niSuS2) + 0{6^). 

In the specific case of the Airy kernels Kj^^^^x, y) := Jj^^ dA Ai(a;+A) Ai(y + A) 
and Kj^^{x,y) := Ai(x + y), 

diKA.isi, S2) + d2KA,{si, S2) = -Ai(si)Ai(s2), (2.22) 

and 

diKA.isi, S2) + d2KA,isi, S2) = 2Ai'(si + S2). (2.23) 
2.2 On the moments (Result [4j) 

Another consequence of the constant shift by a is that all finite moments 
of Ft converge as fast as Without the shift, the first moment would 

converge only as fast as t~^^^, while the variance, skewness, kurtosis would 
of course not be affected by the shift. 

Lemma 2.4. Assume thaE Ft{s) = F{s + f ) + 0{6^)Gt{s) for s e h 
such that F has finite mth moment (with F" G fl C^) and Gt satisfy- 
ing j^(\s\s\^\Gt{s)\ < 00 uniformly in t. Then, 

I s'^dFt{s)= I s™dF(s) + 0(^2) (2.24) 
Jk. Jr 

for all m G N. 

Proof. Let us set It = U /^^ where if" = fl M-t, w = sup /^^ + 5t = inf I^, 
and := + 5t/2. Then, for any m > 1, 

f s"' dFtis) ^6tJ2 ^'><(^) = E ^"(^*(^) - - ^*)) 

= E ^"^(^*(^) - 1) - E + ^triFtis) - 1) + E ^'"^*(^) 

se/+ sei+ s&i- 

-Y.^s + 5trFt{s) + w'^Ftiw - 5t) - w'^iFtiw - 5t) - 1) 

s&I- 

= w-+J2i^--i^ + StDiFtis) - 1) + E (^" - + StDFtis). 

S6/+ s&I- 

(2.25) 



^Note that this condition is stronger than (12.191) and in general not so easy to obtain 
rigorously. 
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Now we set s = s + y so that 



[s + StD = (5 - I)'" - (s + 1)"* = -mStS""-^ + 0{6'f). (2.26) 



Then using Ft{s - |) = + 0{6^)Gt{s) we obtain 
([225]) = w"^ - m5t ^ - 1) - mSt ^ + ^(^^ 



= w"^ - m5t / ds - 1) - m5t / ds + ^(^j') 

^ t« J — oo 

= / s-dF(.) + 0(5n. 
Jr 

(2.27) 

where we used Lemma lA.ll to approximate the sums by the integrals. □ 
2.3 How to fit the experimental data 

For completeness, we explain shortly how to fit the experimental data. We 
partially follow the description of the Supplementary Notes of and use 
their notations. Let us assume that we observed a growth process which is 
thought to belong to the KPZ class. Let S* = e Z, e > 0, be a discrete subset 
of M, where the values of the height function at time t, denoted by ht, lives. 
Let N ':$> 1 the number of experimental measurements and denote by ( ■ ) the 
empirical average over the N experiments. Having f ll.5p in mind, we expect 
to have 

ht^Voot + {Tty/^C + a (2.28) 

where C is a GUE (resp. GOE) Tracy- Widom distributed random variable 
for curved (resp. flat) limit shape, v^o the asymptotic growth velocity and a 
a constant. 

(1) Determine the asymptotic growth velocity Voo- Using 

"^^^'^ ^v^ + ht~^'\ b = ri/3E(C)/3. (2.29) 



dt 

one obtains v^o from the plot (t 



-2/3 d(fet) - 



' dt 



(2) Verify the fluctuation scaling exponent and the fluctuation amplitude 
r. With a log-log plot we can verify if the power 2/3 in 

((/^,-(/^,))2)^(n)2/3Var(C) (2.30) 

holds and at the same time measure the constant T ^ 0. 
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(3) Determine the shift parameter a. Consider the standard KPZ scahng 
ht,resc ■= {ht — Voot)/{T t)^^^. Then, a is measured according to the 
relation 

(Vcsc)-E(C)^a(n)-i/3. (2.31) 

Now that we have determined v^o, F and a, we fit the data vs. the theoretical 
predictions. We set /it,resc := ^t.rcsc — « (Ft)"^/^. 

(4.1) Density: We do a plot of the frequencies of /i^^rcsc with the set 
{S — Voot — ci)/ (Ft)-*^/^ in the abscissa axis. Then we compare this with 
the graph of the Tracy- Widom densites s (-)■ F\s). 

(4.2) Distribution function: We do a plot of the cumulated frequencies of 
ht,rcsc with the set {S — Voot — a)/{Tty^^ in the abscissa axis. Then we 
compare this with the graph of the (shifted) Tracy- Widom distribution 
function s ^ F{s + ^e/{Tty/^)). 



3 PNG and TASEP 

In this section we determine the value of the order 1 shifts for the PNG and 
TASEP models, both with fiat and curved geometry. 



3.1 Flat PNG 

In [H] the formula for the height function ht at time t for the fiat PNG was 
obtaineclfl. It is shown that 

P{ht{0) <H) = det(l - At™),2(|H+i,/f+2,...}) (3.1) 

with 

i^™(xi,x,) = J.,+.,(4t) = ^£dz^,:^, (3.2) 

where J„ is the standard Bessel function (we use the conventions of [I])0. 
As t — )■ cxo, we consider the scaling 

H{s) = 2t + s{2tf''^ en =^ s e It = {N - 2t){2t)-^/^ (3.3) 



^For the one-point distribution there exists also a formulation in terms of Fredholm 
Pfaffian. 

^"With the notation Fs, with S a set, we mean any simple counterclockwise oriented 



path encircling the set S. 
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Under this scaling it is known that [TT] 



^Ai{si + S2)=Kj^,{si,S2) (3.4) 

as t — 7- oo and uniformly for Si, S2 in bounded sets. Moreover, there are ex- 
ponential bounds for the decay of Kf^^l™^ (see, e.g.. Appendix A. 2 of [T3] ) 
which ensures that we can take the limit t — )■ 00 inside the Fredholm deter- 
minant, leading tc("1 

lim PihtiO) <2t + s{2ty/^) = det(l - KA,)Lms,oo)) = Fgoe{2s), (3.5) 

t—^co 

where Fgoe is the GOE Tracy- Widom distribution function [38] . 

Here we focus on the first order correction of K^'^lf^'^ with respect to 
K^-^ and show that it is zero. Since the asymptotic analysis is quite standard 
(see e.g. Lemma 6.1 of [S] for the explanation of general strategy), here and 
in the next sections we indicate only the important steps. 

Proposition 3.1. Uniformly for Si,S2 in a bounded subset of It, 

<r1™^(^l,^2) =/a(^l,^2) + 0(t-^/^). (3.6) 

Proof. We have 

(9fV/^ r p2t{z-z-^) 

K;;r%u s,) -^^i Az (3,7) 



i.resc V-i.-^y 27ri /ro z4t+(si+S2)(2t)i/3+l 

The function z ^ z — z^^ — 2\nz has a double critical point at Zc = 1. The 
steepest descent path can be taken to be coming into Zc with an angle e~'^^^^, 
leaving with an angle e'^^^^, and completed by a piece of a circle around zero 
with a radius strictly larger than 1. Then, the leading term in the asymptotic 
of Kf^l™'^ comes from a t~ ^/'^-neighborhood of Zc- Setting z = 1 + Z{2t)~ 
and doing the large t expansion of the integrand in (13.71) . one obtains 

/ ^3 



1/3 



'"t,resc 



i^'rer^(^i, ^2) = — / dZ exp ( =- - 

2tTI Jooe— i/3 V 3 



^^In [TT] the result is for joint distributions of the height function at different posi- 
tions. The one-point distribution was also obtained through its relation with symmetrized 



permutations [5]. 
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where we have set for simphcity C := Si + S2. Using the contour integral 
representation of the Airy function flA.lip we have 



<t™^(^i,^2) = Ai(C)+t^^/^ 



Ai-(C) CA"(C) Ai^^HC) 
21/3 2'*/'^ 



(3.9) 



Finally, using the identity Ai"{() = ( Ai{Q of the Airy function one readily 
gets that the square bracket is equal to zero. □ 

3.2 PNG droplet 

The formula for the height function ht at time t for the PNG droplet was 
determined in [27]. Let us fix c G (—1, 1). Then, 

nUct) <H) = det(l - K,7-™^),2(|^+i,^+2,...}) (3.10) 

with 

i^rr™^(^i'^2) = ^ j.,+K2tv^r^)j.,+K2tv/r^). (3.11) 

We consider the case c = (and drop the index c) since the general case is 
simply obtained by replacing t by t\/l — (? . An integral representation of 
the kernel is given by 

(27ri)2 Jp^ /p^ ^ ^2t{w-w-^) z-w' 
As t — 7- oo, we consider the scaling 

if (s) = 2t + + a G N ^ s G It = (N - 2t - a)t-^/^ (3.13) 



™«(x„x,) = ^£d«;£ dz "^^^^ (3.12) 



for a t-independent constant a to be specified later. In [27j it is proven that 
the rescaled kernel converges to the Airy kernel, namely 

^™er™^(^i,^2) := t^/=^i^r^™«(^(^i),^(^2)) 

-> / dAAi(si + A)Ai(s2 + A) =/s:^2(si,S2), (3.14) 

as t — oo and uniformly for si,S2 in bounded sets. Moreover, exponential 
bounds for the decay of Kff^^"^^ ensure tha 



lim P(/it(0) < 2t + st'l'' + a) = det(l - i^^jL2((,,^)) = Fgue(s), (3.15) 



^^Thc extension to joint distributions was obtained in [27j . while the one-point result 
is reported in |26j using a mapping to the Poissonized longest increasing subsequence 



problem, which was already solved in [3]. 
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where -Fque is the GUE Tracy- Widom distribution function 

The first order correction of -f^™es™^ with respect to K^^ is the following. 

Proposition 3.2. Uniformly for Si,S2 in a hounded subset of It, with the 
choice a = 1/2, 

^t,rc™^(^i, ^2) = K^M, S2) + 0{t~'/% (3.16) 
Proof. The rescaled kernel is 



+1/3 r r t{z-z-^) 2t+S2t^/^+a-l 1 

^[rr^(5l, ^2) = 7^ i dwi dz 4 Tn — • 

«'=-esc V 1, 2) 27ri 2 / 7 (,t{w-w-^) 2t+siiV3+a 



(3.17) 

Here, we have to integrate over two contours, the one in the z-variable en- 
closing the contour in the u;- variable. The steepest descent path for z can be 
taken as in the flat case, the one for w leaves the critical points with an angle 
g27ri/3 then is completed by a piece of a circle of radius strictly smaller 
than 1. Doing the change of variables 

z = l + t-^/^Z, w = l + t~^/^W, (3.18) 

we eventually get 



i^--™G(si, s,) = I ml dz 



-1/3 



(27ri)2 y^,-2.i/3 y^,-w3 e^'l^-^^^ Z-W 
siZ^ - S2W^ Z^ - 



aZ + {a-l)W 



+ o{t-^'^) 

(3.19) 



where the integration paths do not intersect. At this point we see that setting 
a = 1/2 the first order term is antisymmetric. In particular we can choose 
the paths to satisfy Re Z > Re ly and use -^h^ = dXe~^^^~^^ to get 

^™:sc™''(^i,^2) = K^,is„S2)+t~'/'{Pis„S2)-Pis2,s^))+Oit-'/') (3.20) 
with P given by 



P{s,,S2) = - I dAAi(s2 + A) 



dsi ^dsl dsf 



d^ 

Ai(si + A). (3.21) 



Using Ai"(x) = x Ai(x) and integration by parts one then shows P{si, S2) = 

P(S2,S1). □ 
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Without the shift by a in the scahng, the result would have been 

^t,rc™^(^i,^2) = K^,{suS2) - Ai(sO Ai(s2))) + Oit-'/') (3.22) 

from which we can read off the shift a = 1/2, compare with Remark 12.31 

Remark 3.3. The shift by a = 1/2 is actually independent of c G (—1, 1), 
which is due to the fact that it is built up during the first stages of the growth 
process and for large t it converges to 1/2. Therefore for large t, the shift at 
time ty/l — <? is the same as for the model at time t. 

3.3 TASEP with alternating initial condition 

Now consider TASEP with alternating initial condition, Xk{t = 0) = —2k, 
/c G Z. The joint distribution of particle positions for this initial condition 
has been determined in [TD]. For one particle (here — x„(t) plays the role of 
ht), we have 

P{xn{t) >X) = det(l - <:™P),2(|„„^_2,x-i}) (3.23) 

with 

-^£d^e"-"' (iJ,).„.„.r (3-24) 
As t — )■ oo, we consider the scaling 

X{s) = -In + |t - st^l^ -aeZ s G = (Z - t/2 + a)t-^^^ (3.25) 

for a t-independent constant a to be specified later. It is knowiJ^ that [10] 

<re™''(si,S2) := 2^(^^)-^(^^)t^/3^ftTASEP^^^^^)^^(^^)) 

->JU(S1,S2) (3.26) 

as t — 7- oo and uniformly for Si,S2 in bounded sets. Moreover, exponential 
bounds for the decay of Kffl]'^^^^^ ensure that 

lim P{xo{t) > t/2 - st^/'' -a) = det(l - /U)l2((s,oo)) = Fgoe(2s). (3.27) 

The first order correction of Kf^^JJ^^^^ with respect to Kj^-^ is given as 
follows. 

^•^The prcfactor 2"^(''2)~'^(^i) is just a conjugation, which does not change the underlying 
determinantal point process, but it is needed to have a well-defined limit. 
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-2.0 -1.5 -1.0 -0.5 0.5 1.0 

S 

Figure 1: Distribution function of Xn=[t/4:]{t) for TASEP with alternating 
initial conditions and t = 1000. The number of runs is 10^. The dots are the 
plot of (s G /t,P(x[j/4](t) > X{s))) with 6t = t'^'^ and a = 1/2. The solid 
line is (s, Fgoe(2s + ^St)) while the dashed line is (s, Fgoe(2s + ^6t — aSt))), 
where the shift by a6t for dashed line follows from the definition It, see (13.251) . 




o- 



-2.0 -1.5 -1.0 -0.5 0.5 1.0 

S 

Figure 2: Probabilities for Xn=[t/4]{t) for TASEP with alternating initial con- 
ditions and t = 1000. The number of runs is 10^. The dots are the plot of 
(s e It, P{x[t/4]it) = X{s))6t^) with 5t = t-^'^ and a = 1/2. The solid line is 
(s, 2F{(2s)) while the dashed line is (s, 2F{(2s - a5t)). 
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Proposition 3.4. Uniformly for Si, S2 in a bounded subset of It, with the 
choice a = 1/2, it holds 

Kl:lIJ''^''{Sl, S2) = K^M, S2) + r^/^Kasym(Sl, S2) + Oit^'/'), (3.28) 

where K^symisi, S2) = ^(^2 - sj) Ai(si + S2). 
Proof. The rescaled kernel reads 

^fres?'''''(^i,52) = - — dze'^'-^^ \ \, (3.29) 

The function z t— ?■ 1 — 2^; + ^ In has a double critical point at = 1/2. We 
choose as steep descent path the one coming into Zc with angle e™^'^^ leaving 
with angle e~^'^/^, and continued by a piece of a circle around 1 with radius 
1/2. Setting Z = t'^/'\2z - 1), we get 



iir/3 



ZTTl 



dZ e^'' 



X (^1 + r^/^^ (^(1 - 2a)Z + ^^y^^') + 0(r2/3)^ . (3.30) 

Thus we see that in order to make the first order correction antisymmetric 
we need to choose a = 1/2. With this choice, 

KfS.^^^''{si.S2) = Ai(si+S2) + i(^2-^i) Ai"(.i+S2)r^/=^+0(r2/3) (3.31) 
and the statement follows using Ai"(si + S2) = {si + S2) Ai(si + S2). □ 

3.4 TASEP with step initial condition 

Now consider TASEP with step initial condition, Xfc(O) = —k, k = 1,2, . . . 
The joint distribution of particle positions for this initial condition can be 
found for example (as special case) in [H]. For one particle, we havj^ 

P{xn{t) >X) = det(l - K|7''^'''^).2({.,.,^_2,x-i}) (3.32) 

with 



(27ri)2 Jp^ e^'^yi-wj z''+''^+^ z - w 

(3.33) 



^"^The formula for the one-point distribution can be also given by a last passage percola- 
tion model, which can be analyzed by determinantal line ensembles leading to the Laguerre 
kernel [I^. Joint distributions for the related last passage model can be determined via 
Schur process [HHOlllS] . 
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Now consider a particle that at time t is in the "rarefaction fan", i.e., it is in 
the region with decreasing density strictly between and 1. Such particles 
have a particle number n = at for some a G (0, 1). Thus, we define for a 
couple (n, t) the value of a := n/t and we assume that this value for large 
n and t is clearly away both from and 1. Then, the scaling for t — oo is 
given by 

n = at G N, X{s) = -n + (1 - - sc2t^'^ - a G Z, (3.34) 

so that 

It = {'L-{l-^/aft + a)c2H-^'^ (3.35) 

with C2 = cr~^/^(l — v^)^''^, and for a t-independent constant a to be specified 
later. It is also known that [S] 

j^^stepTASEP / 11/3,-1 /— nX(si)-X(s2) r^stcpTASEP / t;^/ n ( ^ \\ 

-^KaM.S2) (3.36) 

as t — !■ 00 and uniformly for Si,S2 in bounded sets. Moreover, exponential 
bounds for the decay of K^^^^^^^^^ ensure that 

lim P(x„(t) > X{s)) = det(l - KA,)Lms,oo)) = i^GUE(s). (3.37) 

Now let us focus on the first order correction. 

Proposition 3.5. Uniformly for Si,S2 in a bounded set, with the choice 
a = 1/2, it holds 

Kt^:ilsf''''isu S2) = K^,{s,, S2) + c^H-'/'K^,^{s,, S2) + Oit-'/'), (3.38) 
where A'asym(si, 82) = P{si, 82) - P(s2, 81), with 

P(si,S2) = 11 dA Ai(si + A) 
^ Jr+ 

X (^Ai'{82 + A) + S2 Ai"(s2 + A) - ^~jy^ Ai(^)(g2 + A)j . (3.39) 

Proof. With Ci = l — 2y/a we can write X^s) = Cit — C28t^^^ — a. The rescaled 
kernel then reads 

^stepTASEP , , _ C2tV^(l - ^Y.is^-s.)t^^^ 



dz (h dw—{ ^ rj. . (3.40) 
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Figure 3: Distribution function of Xn=[t/i\it) for TASEP with step initial 
conditions and t = 1000. The number of runs is 10^. The dots are the plot of 
(s e Ji,P(x[t/4](t) > X{s))) with 6t = (t/2)-i/3 and a = 1/2. The solid line 
is (s, Fgue(s + k^t)) while the dashed line is (s, -Fgue(s + k^t — a^t)))- 



0.4 



0.3 



0.2 - 



0.1 - 



0.0 




Figure 4: Probabilities for X[t/4]{t) for TASEP with step initial conditions 
and t = 1000. The number of runs is 10^. The dots are the plot of 
(s G Jt,P(x[t/4](t) = X{s))5r^) with 6t = (t/2)-i/3 ^nd a = 1/2. The solid 
line is (3,^2(5)) while the dashed line is (5,^3(5 — a6t)). 
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The function z y-^ z + a\n{l — 2;) — (a + ci) Inz has a double critical point at 
= 1 — ^/a. The steepest descent path for z can be taken such that it comes 
into ^ with an angle e~'^'^^/^, leaves with an angle e^'^'/^, and is completed by 
a piece of a circle around zero of radius strictly larger than ^. The steepest 
descent path for w comes into ^ with an angle e'^^''^, leaves it with an angle 
e"'^^/^, and is completed by a piece of a circle around zero of radius strictly 
larger than 1. By the change of variables 

z = i + c^H-^I^Z, w = i + c^H^^/^W (3.41) 

with C3 = cr~^/^(l — ^/a)'^/^ and a large t expansion of the integrand, we 
have 



-1 pooe ' pooe ' 

<Sf"'(^i'^2) = 7^/ dZ dW 
xfl + C2 4-1/3 f (a- l)Z-aP^+ ~ +c,{Z'-W') ] +0{t-^"'' 



(3.42) 

with C4 = (1 — 2a/o")/(4a/o"). The choice a = 1/2 makes the first order 
correction of the kernel antisymmetric. Finally, we can choose the paths 
satisfying ReZ <ReW and use = dX e"^(^-^) to obtain (IX^ . □ 

Remark 3.6. Looking at Figures [THU one has the impression that TASEP 
with alternating initial conditions is already "closer" than with step initial 
conditions to its asymptotics at time t = 1000, which is confirmed by the data 
in Table [U This is to remind the reader that although in both cases the error 
is (9(t~^/3), depending on the prefactor one still might see some differences of 
the accuracy for not too large times t. The slower convergence for curved vs. 
flat geometry holds also for the PNG model as verified numerically by Richter 
in his diploma thesis [29] adapting the numerical approach of Borneman [7]. 



4 PASEP 

Consider the partially asymmetric simple exclusion process on Z in continu- 
ous time with step initial condition. A formula for the one-point distribution 
of the nth particle from the right has been derived in P^HD] . The expression 
is this time not just a Fredholm determinant, but an integral in the complex 
plane of a Fredholm determinant. A rigorous large time asymptotic analysis 
is in [H], in which it is shown that particles in the rarefaction fan fluctuate 
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TASEP, t = 1000 


Mean 


Variance 


Skewness 


Kurtosis 


Alternating IC 


-0.60495 


0.4027 


0.282 


0.143 


s ^ i^GOE(2s) 


-0.60327 


0.4019(5) 


0.293 


0.165 


Relative error 


0.28 % 


0.18 % 


-3.6 % 


-13 % 


Step IC 


-1.7949(7) 


0.842(2) 


0.19(0) 


0.06(7) 


s Fgue(s) 


-1.77109 


0.8132 


0.224 


0.094 


Relative error 


1.3 % 


3.6 % 


-15 % 


-29 % 



Table 1: Comparison between alternating and step initial conditions. The 
data comes from the simulation used for the previous figures. 



asymptotically according to the CUE Tracy- Widom distribution Fgue- The 
scaling limit to be considered ij^ 

n = ateZ, X{s) = ci{a)t-sc2{a)t^/^ -aeZ (4.1) 

where 

ci(a) = 1-2V^, C2{a) = a~'/%l - (4.2) 

The t-independent constant a will be specified later. Then in [H] it is proven 
that 

lim P(x„(t/7) > A:(s)) = Fgue(s) with 7 = p - g > 0. (4.3) 

t—^oo 

Our result on the first order correction is the following. 
Proposition 4.1. Let p G (|, 1], g = 1 — p, and set 

°° 11 
jr^p -q 2 Vcr ^ 

Then for large time t it holds 

P (a;„(t/7) > X{s)) = Fgve{s + l6t){l + ^(r^/^)), (4.5) 

for s in a bounded subset of It = (Z — Ci{a)t + a)5t with 5t = C2{o')^^t^^^^ . 

Remark 4.2. Note that the previously discussed TASEP with step initial 
conditions is a special case of this result, with p = 1 — q = 1, since ai^o = 



-'^^ There is a minor difference with respect to the papers of Tracy and Widom. To get 
their framework we need to apply the transformation x —x. 
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From this result one can easily get the corresponding result for the height 
function. Let rj^lt) be 1 if there is a particle at site x at time t, and zero 
otherwise. Then the height function is defined by 

r 2J(0,t) + - 2^2/ W)' ^ ^ 1' 

h{x,t)=l 2J(0,t), forx = 0, (4.6) 

[ 2J(0,t) - Ei(l - 2r/,(t)), for x < -1, 

where J(0,t) = YliyyoVyi'^) ■ To get the result for the h from Proposition 14.11 
one simply uses the identity 

P(xj(t) >x) = P(/i(x, t) > 2] + x) (4.7) 

with the following result. 

Corollary 4.3. Let x = E Z and set 

H{s) = 2 (1 + ^')^ - ^ + ^ = 2''^''^ " ^^-^^ 

Then, for large t it holds 

P(/i(x,t/7) > His)) = Fgue(s + i5f)(l + 0(t-'/=^)), (4.9) 

for s G where is defined by the requirement H{s) E 2Z if x is even 
and H{s) E 2^ + 1 if x is odd. Since the lattice width of h is 2, we have 

2^/3(1 -^2)-2/3^-l/3^ 

From this result it follows that the critical value Pc of the asymmetry in 
PASEP such that the density (and the moments) of the rescaled integrated 
current are correct up to order 0{t~'^^^) is the solution of 

ap^,i_p^ = i = 0.7822787862. . . (4.10) 

In Figure [5] we plot the function 2ap.i_p — 1. 

Proof of Corollary \4.!^ Let us define a linearization of the distribution func- 
tions by 



Ftis) 
and similarly 



fp(x„(t/7) > ci(a)t - sc2(a)ti/3 - i + ^ap,^), if s e 
[linear interpolation, otherwise 



(4.11) 



[linear interpolation, otherwise. 
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Then, Proposition 14.11 tell us that 

Ft{s - \bt) = Fgue(s)(1 + 0(r2/3)), s e M, (4.13) 
and we want to show that 

Ft%s - i5f) = Fgue(5)(1 + 0(r2/3)), for s e + \5l (4.14) 
For s G /f + using (14.71) we have 



2^t J ^ V'-'VS'')''/ // — V" 2' 

P(/i(et,t/7) > i^(s) + 1) = F{x^t{th) > ^t) 



(4.15) 



with cr = (iJ(s) + 1 — .^t) / (2t). With this value of a, an algebraic computation 
gives 

= ci{a)t - sc2(a)t^/^ - (a - 1/2) + 0(ri/=^). (4.16) 
Since sc2(cr)t^/3 + (a - 1/2) = (s - |5t)c2(o-)ti/3 + a we get that 

^'(^-I^f) = dllSl) = F,{s-l6t+0{r'/')) = FaMs){l+0{r'^')) (4.17) 

where in the last step we used (14.131) coming from Proposition 14.11 □ 

Remark 4.4. As p — 1/2 the model becomes close to the WASEP studied 
in [2l|3Tl|32l[3l] . In particular, our result matches the limit behavior of [3T] . 
Indeed, when the asymmetry /3 := 2p — 1 — )■ 0, the shift for the fitting of the 
density behaves as 

+ j + (4.18) 

with 7e = —dx ln(r(x)) |_^__^ = —0.57721 56649 . . . the Euler constant, so that 
for the height function the shift is then a = (3~^{\n{2/3) — 7e) — ^ + 



Proof of Proposition 4-1 As for PNG and TASEP, we indicate the main 
steps of the asymptotic analysis to get (14. 5 p for s in a bounded set, but we 
will derive bounds for |s| — )■ oo needed to determine moment convergence. 
Set u = cit — C2st^/^ — a and r = q/p < 1. As shown in [12] . 



P(xn(t/7) >u) = ^ /^(/i;r)oodet(l + /iJ, 



(4.19) 
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0.5 0.6 0.7 0.8 0.9 1.0 

Figure 5: The function p i— 2ap^i_p — 1 for p G (1/2, 1]. It crosses the ordinate 
axis as p = pc = 0.7822787862.'. .. 



where (//; r)oo is the g-Pochhammer symbol (see Appendix [B] for identities) 
and the integral is taken over a circle around the origin with radius in the 
interval (0,r). The operator has kernel 

1 r exp(-^)(l -C)T" fil^A) 



exp(^)(l-V)"(r/0"+^ C-V 



where rj,r]' are on a circle around with radius r G (r, 1) and C runs on a 
circle around with radius in (l,r/r). For 1 < |z| < r~^, the function / is 
given by 



f{f,,z)= J2 1 A-^'' (4-21) 



which extends analytically to C* \ {r^ : k G Z}. 

The function ( i— )■ yi:^ + crlnC + ci ln(l — Q has a double critical point 

at ,^ = — Y3^, and the steepest descent path can be taken such that the 
?7-contour for is a pair of rays from ^ in the directions ±7r/3 completed 
by a circle around zero of radius strictly smaller than 1, and the (^-contour 
is a pair of rays from ^ — in the directions ±27r/3 completed by a circle 
around zero of radius strictly larger than 1. We then do the transformations 

C = ^ + c^H~'/h, 7] = ^ + CsH-'/'w, 7]' = ^ + c^H~'/'w, (4.22) 

with C3 = a~^^^(l — ^/aY^^. Expanding / around z = 1 yields 

/i /(/X, z) = -^ + g{^i) + 0{\z-l\), (4.23) 
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where 

OO ^ OO 

k=0 ^ k=l ^ 

This expansion is obtained by dividing the series into {0, 1, . . .} and 
{. . . , -2, -1}, using 

T^ = r^-1' (4.25) 
1 — r'^jjL 1 — r'^yU 

and then change the variable k ^ —k m one of the sum. After a large t 
expansion, the kernel /xJ^ can be written as 

/.e^WSoo g«}3/3-«i(s+<5t/2) ^{w-z)St/2 

dz ■ 



27ri ye-2.i/3oo e2V3-2(s+<5i/2) (^^ _ y;)^^ _ 

exp((c4^4 - \sw^ - + a)^)c2 H-^/' + 0{t-^/'^)) 



(4.26) 



exp((c4z4 - lsz2 - + a)z)c2 4-1/3 + C(t-2/3)) 
with C4 = ^4^^ . Since Re(2: — w) < 0, we have 

/OO 
dxe("-'")", (4.27) 
-+<5t/2 



3(z-«i)(s+5t/2) 



W — 2; 

and plugging fl4.27p into 04.261) gives 



dx / dz 



27ri y,+5,/2 Je~^-^/^oo e^y^-'^ z-w 

X — ? ^—^ ^ (4.28) 

exp((c424 - |S22 - (M^ - 1 + a)2)c2 W3 + C(t"2/3)) 

The operator is the product Ci C2C3, where the factors have kernels 

1 e-^'/3 



Ci{w,z) 



2m z — 

1 -c^-li-i/S [c424-s22/2-2(3W/v^-l/2+a)] +0(f-2/3) (4.29) 



6*2(2;, x) = —e^^^e 
zvri 

The operator C3C1C2, which has the same Fredholm determinant, acts on 
L2(s + St/2, 00) and has kernel with (x, y) entry given by 

— WX 1 , 

— ^ / dW dw; ^ f 1 + 0(t-2/3) 

(27ri)2 y,-2.i/3^ ie-i/Soo e-V3-.y 2; - w; V 

+ ((c4(w^ - z^) - Isiw^ -z^)-{'-^-l + a){w- z))c^H~'/'\ (4.30) 
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Using again Re{z — w) < 0, we can write = dXe ^\ Thus, ( l4.3Up 
equals 

- K^,{x,y) + c^'r'/\K,,yUx,y) + K,,Ux,y)) + ^(r^/^), (4.31) 
where i^asym(a;, y) = P{x, y) - P{y, x) with 



P{x,y) 



POO 

/ dAAi(a; + A) 
Jo 



Ai(y + A) 



(4.32) 



is asymmetric, and Ksyra{x,y) = — — | + a) Ai(x) Ai(?/) is symmetric. 
Hence, we have 



1/3 



+ 0(r2/3)). (4.33) 



with Xs the projection onto (s + (5(/2, oo) and 

G = 7^J -^?(/i)(/i;r) 



(4.34) 



We will show that G = ttp^g so that by choosing a = j — -^cip,q the prefactor 
of the first order correction vanishes. First note that 

-i./ ^(/.;r)^f;-^ = 0, (4.35) 
as the integrand has no poles inside the unit circle. So, we have 



oo 



We use 



1 



1 1 
- + 



(4.36) 



(4.37) 



the fact that (yu; r)oo is analytic inside the integration domain, so that the 
sum of the contributions of the simple poles gives 



(4.38) 



k=l 
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There is a simpler expression for G. Using the identity flB.2p we get 



G 



EE 

k=i e=i 



-1)V 



_ki 



E 



£=1 



-l)V 



e e{e~i)/2 e 



r;r 



(4.39) 



Then we use ( IB.3P and (a; r)f/(ar; r)^ = {a — l)/(ar^ — 1) to get 



r;r - 1 



hm — ^( i0i( " 



hm — ^ hm \2<Pi\^'^ 

a-5>l a — 1 ^^oo V V <^''" 



(4.40) 



where we used ( 1B.4I) in the last equality. Finally, the g-Gauss identity ( IB.SP 
leads to 



G = iim^iim 



^^(ariT)^ 



a-^i a - 1 /3-^oo V(ar; r)oo(r//3; r)oo 

oo 

5c,ln[(ar;r)oo]|^^-^ = -(9a^ln(l - ar 



r; r 



^=0 



1 ' /oo 

oo 



5:t 



Replacing t = q/p leads to G = a^^g. This and fl4.33p shows that 

P(Xn(t/7) > = i^GUE(s)(l + 0(r2/2)). 



(4.41) 

(4.42) 
□ 



A Discrete sums versus integrals 

Lemma A.l. Let / : M" — )■ M &e a smooth function with djdkf G L^iW^), 
for all j, k = 1, . . . ,n. Then, for d > (small) 

fi^S)- [ d-xfix) =0{6')f2 I d"x|W(x)|. (A.l) 



Proof. We first rewrite 



f d"x/(x)= V / d" 



y/(x5 + y) (A.2) 
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and use Taylor development 

f(xS + y) = f{xS) + djf{xS)yj + Viy^RjA^^ v)^ (^-3) 

with 



2 

3=1 j,k=l 



\R,,k{S,y)\< max \d,dkfix5 + u)\, for y G [-|, f]'^ (A.4) 
to obtain (after integrating over y), 

(A.5) 

The statement then follows because the sum over x converges, as 5 —t- 0, to 

E-fc=i/M.d"^|5A/(^)l- □ 

Lemma A. 2. Let f{xi, . . . ,Xn) = det{K{xi,Xj))i<ij<n with the kernel K 
satisfying 

max{|ir(xi,X2)|,|Wxi,X2)|,|a,a,i^(a;i,X2)|} < Ce-^(-^+^^) (A.6) 
for all Xi,X2 G (s, oo), i,j G {1,2} and some positive constants c, C. Then, 

n 

fc=i 

for all 1 < i, j < n. 

Proof. Let Ki^ (resp. Kj) be the matrix Xj))i<jj<„, with the ith row 

(resp. the jth column) replaced by its derivative w.r.t. the first (resp. the 
second) variable. Then, 

difixi, ...,Xn) = det Ki^ + det K^i (A.8) 

and from this 

didjf{xi, . . . , Xn) = det Kij^ + det Kij + det Kj^i + det Kji (A. 9) 

with Kij^ = {Ki)j^. By Hadamard's bound, the absolute value of an n x n 
determinant with entries in the closed unit disk is bounded by n"/^. It then 
follows 

n 
k=l 

for any 1 < j < n. □ 
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Lemma A. 3. For the Airy kernels, Kj^^^i^^ u) = /m Ai(x + A) Ai(?/ + A) 
and Kjs,^{x^y) = Ai(x + y), assumption l(A.6\} of Lemma \A.^ is satisfied. 

Proof. It is easy to see from the integral representation 

Ai(x) = — / d2e"'/='-^" = — / d2e^'/=^-^", £>0, (A.ll) 
27ri Jooe— i/3 27ri 



that for any 5 > 0, there exists a constant G (0, oo) so that 

max{|Ai(x)|, |Ai'(x)|, |Ai"(x)|} < C^e"^" (A.12) 

uniformly in x G M. 

In the case K{xi,X2) = J^^ dX Ai(xi + A) Ai(x2 + A) we get from the 
bounds flA.12p with e = \, after integration with respect to A, that 

max{|i^(xi,X2)|, \d,K{xuX2)l xs)!} < Ce-^^^+^^^Z^ (^.13) 

for some constant C > and all i,j G {1,2}. 

The case K{xi,X2) = Ai(xi + X2) is even easier, since the bound ( IAT3|) 



comes directly from ( 1A.12|) . □ 

B g-Pochhammer symbols, g-hypergeometric 
functions 

Here we collect some identities on g-Pochhammer symbols and g-hyper- 
geometric functions, used for the PASEP. We use the standards as in [23]. 
The g-Pochhammer symbol is defined by 

00 n— 1 

{^^]q)oo = \{{l- m^). and {^l]q)n = \{{l- m^)- (b.i) 

fc=0 fc=0 

They satisfies the following identities: 

(ii-n\ °° 'l'^n"("-l)/2 

1/^9"; <i)oo (g; g)n 

so that in particular (0; g)oo = 1 and (1; g)oo = 0. 
The g-hypergeometric function is defined by 



f ai, ■ 


. . , Ctf 




\bi,. 


■ ■,bs 


g; 



l+s—r 



g)n • • • (Qr; g)n (( _l\n n{n-2)/2y 

(B.3) 
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In particular, it holds 



r-Ws 



Oi, . . . , ttr-l 
bi,...,bs 



lim 



Oi, . . . , 
bi,...,bs 



The g-Gauss identity is 



2V1 



fa, (3 




9; 7^ 


V 7 



(7/«;9)oo(7//5;g) 00 

(7;?)oo(7/(a/3);9) 00 



(B.4) 



(B.5) 
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